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Abstract 

The anticommuting analysis with Grassmann variables is applied to the two- 
dimensional Ising model in statistical mechanics. The discussion inclndes the trans¬ 
formation of the partition fnnction into a Gaussian fermionic integral, the momen- 
tnm-space representation and the spin-fermion correspondence at the level of the 
correlation fnnctions. 


There are many remarkable interrelations and analogues in studies of model systems 
in quantum held theory and statistical mechanics. The two-dimensional (2D) Ising model 
(2DIM) may be a good example of this kind. In its original formulation, this is a discrete- 
spin lattice model of a ferromagnet [1-7]. The important role of 2DIM in the theory of 
phase transitions is due to the possibility of exact solution at all temperatures (yet in a 
zero magnetic held) [1-23]. A remarkable feature is that 2DIM can as well be reformu¬ 
lated in terms of free Majorana helds on a lattice [8-15]. This can be done most ehectively 
in terms of anticommuting Grassmann variables and integrals [8,24]. In this report we 
shortly review the application of anticommuting integrals to 2DIM (also see [18]). The dis¬ 
cussion includes the transformation of the partition function Q into a Gaussian fermionic 
integral, the momentum-space analysis, and the spin-fermion correspondence at the level 
of correlation functions. The mirror-ordered factorization scheme we apply for translating 
Q into Grassmann variable integral was originated in [12-14]. For recent developments in 
the theory of 2D Ising model also see [16-23]. 

In physical aspect, Grassmann variables are in essence the classic (nonquantum) 
fermions [8,9,24]. The rules of integration over Grassmann variables were invented by 
F.A. Berezin in early 60ies in the quantum-held theoretical context [8]. The original 
motivation of his work was in looking for a suitable fermionic counterpart for ordinary 
variables of commuting analysis (bosonic helds) in path integral [8]. Since that time, the 
anticommuting integral calculus has found numerous applications in quantum held theory 
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and condensed matter physics. Before proceeding directly to the Ising model, we shortly 
remember the basic rnles of anticommnting integration. 

In mathematical aspect, Grassmann variables may be viewed simply as formal purely 
anticommuting fermionic symbols [8,9,24]. Given a set of N independent Grassmann 
variables ai, 02 , 03 ,..., qn, one has: OjOj + ajQi = 0 , = 0. The variables may be mul¬ 

tiplied by each other and by complex numbers in a usual way. The linear superpositions 
of Grassmann variables are again purely anticommuting, their squares are zeros. The 
first important identity of the anticommuting analysis follows by multiplying the linearly 
transformed variables: 


N 

6i62&3----&Ar = deti ■ aia2a3...aAr , h = (1) 

Berezin’s rules of the integration for one variable are [ 8 ]: / daj ■ aj = 1, / daj ■ 1 = 
0. In multiple integral, the differentials are again anticommuting with each other and 
the variables [ 8 ]. Any natural function of Grassmann variables oi, 02 ,..., oat may be 
in principle represented as a hnite polynomial in these variables. The integration then 
reduces to the repeating application of the elementary rules [ 8 ]. The Gaussian fermionic 
integral of the hrst kind with integration over complex-variable (Dirac) fermionic helds 
is expressed as the determinant of the matrix from the quadratic fermionic action in the 
exponential: 


N N N 

/ n da*daj exp ( ^ ^ = det A . ( 2 ) 

j=l i=l j=l 

The appearance of the determinant in m may be traced back to (CD- The exponential 
in Q is in fact a finite polynomial due to nilpotent property of fermions. The same 
polynomial follows by multiplying elementary factors exp = 1 -|- ttiAija*. The 

Gaussian fermionic integral of the second kind, for real (Majorana) fermionic helds, is 
related to the Pfafhan: 


N N 


dttN ... da 2 dai exp UEE ttiAijttj'^ = PfaffA, 

2 i=ij=i 


( 3 ) 


where the matrix A is assumed skew-symmetric, A A^ -|- A^A = 0 , in components, Ajj -f 
Aji = 0 , A-- = 0. The number of variables N is assumed to be even. In mathematics, the 
Pfafhan is a combinatorial polynomial associated with a triangular array of elements Ajj . 
In given context, the triangular array is to be identihed with the above-diagonal array of 
the skew-symmetric matrix A from the action. The combinatorics of the Pfafhan is in 
fact identical to that of the fermionic version of Wick’s theorem. 

The 2DIM can be solved exactly only for the homogeneous (translationally invari¬ 
ant) lattice, when all couplings at the congruent bonds are equal. The fermionization 
itself, however, can be performed equally well in the general inhomogeneous case [12,18]. 
Therefore, we start here with a generalized formulation of the model assuming arbitrary 
inhomogeneous distribution of the bond coupling parameters over a rectangular lattice 
net. The Ising spins, Umn = ±1, are disposed at the lattice sites, mn, with m, n = 1,..., L, 


2 


running in horizontal and vertical directions, respectively, where L is the length of the 
lattice side. At final stages N = ^ oo. The hamiltonian is: 


l3H(o)= E E [*: 


d) 

m+ln 


, 7 ( 2 ) 

^mn^m+ln \ ^mn+1 ^mn^mn+1 


m=l n=l 


(4) 


where (3 = 1/kT is the inverse temperature and are dimensionless bond cou¬ 
pling parameters. For finite L, we assume free boundary conditions: cxi+in = = 0. 

The partition function is: Z = H exp where the sum is taken over spin conhgu- 

rations provided by cr^n = ±1 at each site. A typical Boltzmann weight is: exp (baa') = 
cosh 6-|-(ja' sinhfo, since (aa') = ±1, so the problem may be readily reformulated in terms 
of the reduced partition function: 

L L 

Q= Sp {n n( 1+ tm+ln^mn^m+ln) (1 T tmn+1^mn^mn+l) J’ ) (5) 

{^) m=l n=l 


where = tanh b^^^ and Sp(^) stands for the correctly normalized spin averaging, 
such that Sp (1) = 1, Sp (amn) = 0 at each site. 

To convert Q into a Gaussian fermionic integral, the starting point is the fermionic 
factorization of the local Boltzmann weights iniini. We introduce a set of the totally 
anticommuting Grassmann variables, amn, Cimn) ^mn, bmn, ^ bond, and write: 


1 T ^m+ln^mn^m+ln J ® T ®mnC’’mn) (1 T ^m+ln ^mn'^rn+ln) , 

1 T ^mn+l'^mn'^mn+1 J ® A ^mn'^mn) (1 A tmn+1 ^mn^rnn+l) ■ (6) 


In a conventional notation, the bond Boltzmann weights are now presented by doublets 
like and standing under the sign of the diagonal Gaussian averaging 

arising by factorization, where separable factors Amn, Bmn, A^^, B^^ (to be called shortly 
Grassmann factors) are to be identihed from (IH)). At the next stage, we keep in mind 
to group together, over the whole lattice, the four factors with the same spin (the same 
index mn) and to sum over amn = A1 in each group of factors independently thus passing 
to a purely fermionic expression for Q. The Grassmann factors are in general neither 
commuting nor anticommuting with each other, so a special arrangement of factors in 
their global products is needed. The correspondent ordering procedure is called mirror- 
ordered factorization [121 HE] The final result for the factorized density matrix is: 


L L 


L L 


<3(^) = / n n damndamndbmndbmn ^Xp { ^ ^ [ OmnCn A } 


m=l n=l 


m=l n=l 


n 

L L 


Gn[n(i+ b[rm^m-ln'^rnn)(^ A A ttmn'^mn )-n (lA bmn^mn )]}• (7) 


n=l m=l 


m=l 


The averaging over spins is a step-by-step procedure realized at the junction of the two 
m-ordered products in dZD, with given n. At each stage, the averaging over amn = A1 
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of the product of four neighbouring factors ^mn^mn^mn^mn junction results a 

simple polynomial even in Grassmann variables, which can be readily transformed into a 
Gaussian exponential. This corresponds to the local mn-term in the fermionic action in 
the second line of (jHI) below. In this way, one hnally eliminates all spin variables. The 
partition function then appears in the form a Gaussian fermionic integral: 


^ 1 j 1 j 1 j 1 ^ 

Q / 11 11 ^ ^ ^ ^ [ ^mrflmn ^mn^mn 


m=l n=l 


m=l n=l 


1 y-(l)y-(2) ^ * u* I 4-/(2)/j* \(n -\-h ^ -\-n h 11 

' ^mn^mn ^ra—ln^mn—1 ' \^mn^m—ln ' ^mn^mn—1) v^mn ' ^mnj ' ^mn^mn \ J 


( 8 ) 


The above expression is exact, and completely equivalent to assuming free boundary 
conditions for fermions: = 0. The problem is now reformulated as a theory of 

free fermions on a lattice. For further possible transformations of the integral (jHI) and the 
continuum-limit interpretation of the model also see [miiHi. 

For the homogeneous (translationally-invariant) lattice, we put = tl^ 2 ^ for def¬ 
initions see ®. The integral (jHI) can then be explicitly evaluated by passing to the 
momentum space for fermions, which results Onsager’s expression for the free energy per 
site of inhnite lattice (as oo). In the momentum space, the integral (jHj) becomes: 


L-lL-l 


L-lL-l 


<5 = / n n [HY. [%q(^pq + 'bpqbpq + %q^ 


L—pL—q 

p =0 q =0 p =0 q =0 

+ tit2e^ L * L 4 -(ti e* ^ a*q + t2e'' ^ + bpq) | , 


+ 


(9) 


where 6^^, b*q are the new variables of integration introduced by discrete Fourier 

substitution with periodic boundary conditions (we may change boundary conditions into 
most suitable form in view of hnal limit cx) which abolishes boundary effects). 

After a proper symmetrization of the fermionic sum with respect to p,q L — p, L — q, 
the integral decouples into a product of elementary integral factors with four pairs of the 
conjugated fermionic variables per pq site. These integral factors, can be readily evaluated 
by the standard rules of fermionic integration and their product gives the squared partition 
function [niiiHj. The result is: 




n n [(i+^i)(i+^ 2 )- 2 ti(i-ti) cos^ 

p=0 9=0 ^ 


2 * 2(1 



( 10 ) 


which is the exact solution for in the limit L'^ —>■ cx). The correspondent free energy 
per site, -/?/q = logQ | ^ 2 ^^ , then follows: 


277 277 


(^fQ = 


0 0 


dp dq 
271 277 


In 


(1-I-f^)(l-I-^2) — 2fi(l — ^2) cosp — 2^2(1 — cosg . (11) 
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Respectively, the true free energy per site, for Z = E exp(—/3if), can be recalculated 
from Z=( 2 cosli 6 i cosh 62 )^ Q, and we gain: 


-/^fz 


ln 2 + 


27r 27r 



dp dq , 

-In 

2tt 2tt 


0 0 


cosh 2bi cosh 262 — sinh 2bi cosp — sinh 262 cos q 


( 12 ) 


which is the famous Onsager’s result, see Eq. (108) in pp. From the exact solution for 
the free energy, it follows that the point of phase transition in the ferromagnetic case is 
given by the condition 1—ti—^ 2 —ti ^2 = 0, or sinh26i - sinh 262 = 1- The singularity in the 
specihc heat near Tc appears to be logarithmic. For further details on the thermodynamics 
of the two-dimensional Ising model see pm. 

In general, the studies of the correlation functions are important in physical applica¬ 
tions and in the theory of the critical point. Besides, the perturbative expansions produce 
the correlations. The fermionic interpretation of the two-spin correlation function in two- 
dimensional Ising model is also possible but is in fact not so simple since it is nonlocal.^ 
The general statement is that spin-spin correlator {amnO'm'n') is always representable as 
Toeplitz like determinant of large size which elements are all expressible in terms of the 
fermionic Green’s functions like This is in a precise agreement with the results 

of previous combinatorial studies of correlations in 2DIM [21 m (fh® difference is that we 
use anticommuting algebra instead of ordinary combinatorics in due course of derivation). 
The spontaneous magnetization then follows via the Szego-Kac theorem mm 

There are no magnetic field source terms in the fermionic action in (jHl), but at the 
spin level it is possible to express the two-spin function like (ctoctr) in terms of a string 
of local perturbations according to the rule: (ctocth) = (ctocti)(cti(T 2 ) ... (cTi?_iCTR), that 
follows from ct| = 1 , where we assume a local (lexicographic) enumeration of spins by 
travelling from site 0 to site R along some path on a lattice. Now, consider the case of 
two spins disposed on the same horizontal row of a lattice at distance R from each other: 


{(^mn^m+Rn) ( (^mrt^m+ln) (^m-|-ln*^m-|-2n) • • • (^m+R—In^m+ijn)) • 


(13) 


In turn, each entry (Jmn<^m+in may be expressed as a modification of the bond cou¬ 
pling parameters in the partition function since (amnO'm+in)(f + h o'mnO'm+in) = G (1 + 
o'mnO'm+in), where use is made again from (o'mnO'm+in)^ = +1- The problem thus 
reduces, in principle, to the evaluation of the fermionic integral (jH)) with modihed bonds 
along the path. In practical aspect, though, it might be suitable to diminish the number 
of the modihed fermionic terms as much as possible. This may be realized by relating 
o'mnO'm+in dhectly to amnO^mn notice that introducing prefactor amnO-mn under the 

integral in (jH)) effectively annihilates the weight (l+tiamnO'm+in) in the partition function. 
This can be most readily seen introducing amnO^mn under the integral at the level of local 

^ The nonlocal feature in the spin-fermion correspondence at the level of partition function in 2DIM 
is somewhat masked since the Ising spins in 0 and, respectively, fermions in ©, are interacting only 
with nearest neighbours. The nonlocality can also be related more directly to fermionic algebra (anti¬ 
commutativity) if one starts to derive the fermionic representation for the spin-spin correlator from the 
factorized representation for the density matrix 0 . This version is not considered in more detail in 
present discussion. 
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fermionic factorization in This results the correspondence: 


1 T tiCffinCjn+lr, 
The inverse correspondence readily follows: 


1 


( 14 ) 


1 

tl 


1-^1 

tl 


^mn^mn GXp •[ (1 t^') ’ 


(15) 


which enables one to obtain the two-point spin correlation function in terms of a nonlocal 
fermionic string of perturbation under the averaging, to which one can apply fermionic 
version of Wick’s theorem. This eventually results a Toeplitz determinant in terms of the 
fermionic correlations, see m below. It will be suitable yet to express the local terms 
from m as linear fermionic forms. Introducing at each site a pair of auxiliary fermions 
Cmn, Cm we write: 


1 

h 




'^mnO‘rnn 


d,C* d,C -I- (7 ) 


1 * - 1 

_ ^ * -I- _ n 

, ^mn ' , mn 

tl tl 


)}■ ( 16 ) 


By substituting into the string product m the correspondence dl^, one obtains the 
correlator as a nonlocal fermionic string average: 


R-l 


R-l 


{^mn^m-\-Rn') ^ ^rn+an^rn-\-a-\-ln^ ( n [^1 ] + (^1 '^rn+anttm+an)') (l^^) 


Q:=0 


q ;=0 


while with the use of (cni), one gains representation: 

{(^mn(^m+Rn) {cr) ~ i_i^rnn^mn)i^rn+ln^m+ln) ■ ■ ■ i^m+R-ln^m+R-ln)'} 
= det { }i a,f5 = 0,1,.. 


-^mn ^mn "i* ^mn i 




A* = Jl Id_ 

^mn • ^mn ' » 

tl tl 


{^m+an^m+fSn) 7^ 0 ) (^m+on^m+/3n) ~ {^m+an^m+0n) ~ 0 ' 


(18) 


That the correlator in the first line appears to be the determinant (but not the Pfaffian) 
follows from the selection rules outlined in the last line of (HHD. Due to these rules, which in 
turn are based on {amnO^m'n) = (Omn®m'n) = 0 for the correlators on the same horizontal 
row, the Pfaffian recursion of Wick’s theorem reduces to the determinantal recursion.^ 
The determinant that appears in m is Toeplitz determinant of R x R matrix with 


^ The equations {amndm'n) = (Cn®m'n) = 0 follow by a direct calculation via the momentum space, 
cf. (12011 - (12ID . The property that {amnCim'n) = {KnnKn'n) = 0 is a particular feature of a situation under 
the consideration, being valid only for the fermions disposed on the same horizontal row. In a more 
general case, for spins at arbitrary positions, the result will be the Pfaffian (it is known though that the 
squared Pfaffian is again the determinant). 
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matrix element A{a — (3) = (y4m+an^m+/3n)) which in what follows we rather assume 
the notation A{m — m') = {A^^A*^,^. In terms of the original fermions: 


(A 


A* 

mn^m'n 


> = f '5„ 


l-tj 


( 19 ) 


The explicit solution for the only nonzero correlator we need follows directly from 


^ L-lL-l 


2n27T 


PQ PQ 




L^oo 


p=0 q=0 

eipim-m')+iq{n-n') J _ ^2^ ] 


0 0 


dpdq 

(27r)2 (1 + tf)(l + tl) — 2ti (1 — tl) cosp — 2^2 (1 — tf) cosq 


( 20 ) 


which we have to substitute into dm . taking the case n = n'. This results the integral 
for (AmnA^/„), in which the integration over dq can be performed by residues. After 
some straightforward though lengthy transformation in the numerator, and changing p 
2n — p, the hnal expression for the matrix element of the Toeplitz determinant appears 
in the form: 


Stt • • ... 

(A /I* \= f dp -(1 - t2 - A - tit2 (1 + ^2 - A + A^2 -ip(m-m') ( 2 ^) 

\ mn m'n) ^21, | (1 - fs " 6*^) (1 + ^2 " A 6*^ + 6*^) | ^ 


Equivalently, assuming A(m — m') = {AmnAdyrHr^^ the above expression may as well be 
represented in the form: 



2p(m—m^) 


(1 e*P) (1 -fr^^2e“*^) 

(l-fif^e-*P)(l-fr^t| e+*P) 


271 


f ^22) 

J 271 


where the Kramers-Wannier conjugate parameters are: 


t 


* 

1 


i-h 

1 -i- ti 



1_^ ^ p-262 

l + t2 


(23) 


This is a self-dual transformation: t^* = ti, = ^ 2 - The condition for the ferromagnetic 
critical point, 1—ti—^ 2 —^ 1^2 = 0, is equivalent to ti = ^ and ^3 = The matrix 
element A{m — m') depends only on the difference of the indices, for which reason the 
correspondent determinant is called Toeplitz determinant El 13 Notice that A{m — m') 
is yet not a symmetric function with respect to changing the sign of m — m'. The last 
equation in (0^ is merely to dehne the momentum-density function c(p). This function 
c(p) being in essence a phase factor possesses interesting properties like |c(p)p = 1 and 
c{—p) = l/c(p).^ 

^ It is not yet well understood, beyond formal calculation, for what reason one might expect for c(p) 
to possess such particular properties. In any case, the form of c{p) is seemingly related to the special 
role of the zero momentum modes (both ferromagnetic and antiferromagnetic) in path integral 0 with 
respect to the formation of the spontaneous magnetization, M, which feature can in fact be recognized 
in the explicit solution for M itself, see (EU below. For this question also see M 
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The correlator {(Tran^'ra+Rn) is thus expressed in terms of Toeplitz determinant which 
matrix elements A{m — m') are fermionic Green’s functions explicitly given in (j22|l . The 
similar Toeplitz determinants typically appear in classic combinatorial treatments El El 
The analysis (rather complicated) of the Toeplitz determinants of this kind provides the all 
known information, at least, in that part of it that can be derived from the hrst principles, 
about the real-space asymptotic behaviour of the two-point spin-spin correlation function 
in 2DIM El 013 • 

The solution for spontaneous magnetization also follows from the Toeplitz determinant 
j2j. In fact, the squared spontaneous magnetization, M^, is the limiting value of two-spin 
correlator {o'mnO'ra+Rn) as i? —cx), and may thus be obtained by means of the Szego-Kac 
theorem [2] as the limiting value of the Toeplitz determinant resulting in (ITHj) - (E^ as 
i? —> cx). Though the Szego-Kac theorem by itself is rather a complicated mathematical 
statement, as regards its justihcation, see a discussion in [2j, its application to the evalu¬ 
ation of is quite simple. The Szego-Kac expresses the limiting value of the Toeplitz 
determinant in terms of the real-space Fourier components of the function Inc(p), rather 
then c(p), where c{p) is the Fourier image of A{rn — m!) in notation of (I22j) . This change 
from c(p) to Inc(p) simplihes the situation drastically. The spontaneous magnetization 
hnally appears in the form j2]: 


M 




1 


1/8 


sinh^ 2bi sinh^ 2^2 


(24) 


The above expression is valid in the ordered phase, in the domain sinh26i •sinh262 > 1- It 
is interesting that despite of the fact that the hnal expression for the spontaneous magne¬ 
tization (E^ is seemingly quite simple, it is even not the elliptic integral, its derivation by 
any known method is few times more complicated than the evaluation of the free energy 
and specihc heat in the correspondent approach. The fermionic analysis, though simplihes 
combinatorial treatment, demonstrates the same feature. Respectively, the problem of a 
better physical understanding of the ordering phenomena in 2DIM in terms of fermions 
still remains. What is ordered in two-dimensional Ising model in terms of fermions? In 
this respect, it might be worth mentioning that there are two other important problems 
with fermions in two dimensions, the quantum Hall effect and the high-Tc superconduc¬ 
tivity in copper oxides. Finally, the two-dimensional Ising model, within the fermionic 
interpretation, may be considered in a common range with other typical problems in 
quantum statistics and quantum held theory. 
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